A new inverse quasifission mechanism to produce neutron-rich transfermium nuclei 



David J. Kedziora^ and Cedric Simenel^' ^'0 

^Department of Nuclear Physics, Research School of Physics and Engineering, 
Australian National University, Canberra, Australian Capital Territory 0200, Australia 
^CEA, Centre de Saclay, IRFU/Service de Physique Nucleaire, F-91191 Cif-sur-Yvette, France. 

(Dated: March 11, 2010) 

Based on time-dependent Hartree-Fock theory, a new inverse quasifission mechanism is proposed 
to produce neutron-rich transfermium nuclei, in collision of prolate deformed actinides. Calculations 
show that collision of the tip of one nucleus with the side of the other results in a nucleon flux toward 
the latter. The role of nucleon evaporation and impact parameter, as well as the collision time are 
discussed. 



o 



(N 
> 

cn 
o 
o 



I. INTRODUCTION 

The quest for, and study of the heaviest elements has 
involved much experimental and theoretical effort in re- 
cent years. Their existence relies only on stabilizing 
quantum shell effects, which make them ideal to test 
quantum mechanical nuclear structure models. Both the 
location of the predicted island of stability in the su- 
perheavy element (SHE) region |l|-[3l and spectroscopy 
of transfermium nuclei {Z > 100) [1, [Bj] are needed to 
constrain these models. A natural way to form such nu- 
clei is through fusion of heavy nuclei, followed by neu- 
tron and gamma evaporation from the compound nu- 
cleus. SHEs have been produced either in "cold" fusion 
reactions based on closed shell target nuclei [l, ? or in 
"hot" fusion reactions involving actinide targets [H, [9|. 
The heaviest element, containing 118 protons, has been 
synthesized with the latter technique [s]. However, a- 
decay chains of SHEs formed in hot fusion end in a region 
of unknown neutron-rich isotopes with 104-110 protons. 
Thus, it is necessary to study this region of the nuclear 
chart in order to provide a better identification of the 
decay daughters and confirm these SHEs. 

Fusion-evaporation cross-sections decrease rapidly 
with the product of the charges of the reactants, down 
to few picobarns for SHEs. These cross sections are too 
small to allow a detailed study of nuclear structure. For 
instance, a basic property like the mass has been mea- 
sured only recently for 252-254]\Jq fusion products with a 
Penning trap mass spectrometer p^]. Furthermore, fu- 
sion reactions lead usually to neutron deficient compound 
systems, which, in addition, decay by neutron emission. 
It is therefore worth exploring other reaction mechanisms 
to produce and study the heaviest nuclei, and, in partic- 
ular, their neutron-rich isotopes. 

An alternative way to form neutron-rich heavy nu- 
clei is to consider multinucleon transfer in such a way 
that one eject ile gets heavier that any of the reac- 
tants This process is sometimes called "asymmetry- 
exit-channel" [l2[ or "inverse" [l3| quasifission, as the 



mass asymmetry of the outgoing fragments has increased, 
whereas "standard" quasifission tends to reduce this 
asymmetry. Such a process has been investigated exper- 
imentally considering either a lig ht- medium mass pro- 
jectile on an actinide target |l3-[l8j, or actinide colli- 
sions [l9l-[23j. Recent theoretical studies of multinucleon 
transfer have been performed in the dinuclear system 
model (DNS) ^S^lsj? using multidimensional Langevin 
equations [13^, 25-^, in the constrained molecular dy- 
namics model [28], in the improved quantum molecu- 
lar dynamics approach [29|, |30|], and within the time- 
dependent Hartree-Fock (TDHF) theory [31]. In partic- 
ular, it is predicted that shell effects in the ^^^Pb region 
should favor inverse quasifission [l3|, l20|, |2^ . Indeed, as 
one actinide falls into the valley of the potential energy 
surface toward the magic numbers Z = 82 and = 126, 
the mass and charge of its collision partner increases cor- 
respondingly. 
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Multinucleon transfer depends also strongly on defor- 
mation and relative orientation of the nuclei [30|, Isij . 
This should play an important role in actinide collisions 
as nuclei have strong prolate deformations in this region 
of the nuclear chart [32]. In particular, it has been shown, 
in the case of the symmetric central collision ^^^U+^^^U, 
that a nucleon flux appears in the neck when one nucleus 
has its deformation axis aligned with the collision axis 
and perpendicular to the deformation axis of the colli- 
sion partner [31]. In this case, the "aligned" nucleus loses 
nucleons. The main goal of the present paper is to inves- 
tigate a new inverse quasifission mechanism due to such 
orientation effect in initially mass and charge asymmetric 
collisions of actinides. As an illustration, we perform cal- 
culations for the Th+2^°Cf reaction within the TDHF 
framework. In section [III we present briefly the TDHF 
theory and give some numerical details of the calcula- 
tions. Then, the results are presented and discussed in 
section UlTl Finally, we conclude in section HVl 
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II. THE TIME-DEPENDENT HATREE-FOCK 
APPROACH 

A. Theory 

The TDHF theory has been proposed by Dirac 
and apphed in nuclear physics [34, 35], including actinide 
collisions [36|, with Skyrme effective interactions [37]. In 
its Liouville form, the TDHF equation is written 

in^ = [h[p],p]. (1) 

It gives the evolution of the one-body density matrix p 
assuming that the system is always described by an an- 
tisymmetrized independent particle wave function to en- 
sure an exact treatment of the Pauli principle during time 
evolution [38] . The one-body density matrix can be used 
to compute expectation values of any one-body observ- 
able and its evolution, within TDHF, accounts for one- 
body dissipation mechanisms. The latter are known to 
drive low-energy reaction mechanisms as the Pauli block- 
ing prevents nucleon-nucleon collisions. 

The one-body density matrix of an independent parti- 
cle state can be written 

p{rsq,r's'q')= ^ V*(r'sV)^iM, (2) 

where {(pi} are the occupied single particle wave func- 
tions, Ai and A2 are the number of nucleons in each nu- 
cleus, and r, 5, and q denote the nucleon position, spin 
and isospin, respectively. The single-particle Hartree- 
Fock Hamiltonian h[p] is self-consistent and can be ex- 
pressed as 

where E[p] is the Skyrme energy density functional 
(EDF) modeling the interaction between the nucleons. 
The EDF is the only phenomenological ingredient, as it 
has been adjusted to reproduce nuclear structure proper- 
ties [39|. In practice, the TDHF equation ([1]) is written 
as a set of nonlinear Schrodinger-like equations for the 
occupied single-particle wave functions 

in^^ = h[pmy,,{t). (4) 

Realistic TDHF calculations in 3 dimensions are now 
possible with modern Skyrme functionals including spin- 
orbit term (4Q|-|43| and supercomputers allow simulation 
of realistic actinide collisions (sij . 

B. Numerical details 

The nuclei are assumed to be initially in their Hartree- 
Fock (HF) ground state. The HF and TDHF calculations 



are both performed with the SLyM Skyrme EDF [ioj , 
allowing for a fully consistent treatment of nuclear struc- 
ture and dynamics. HF ground states are generated by 
solving the stationary version of Eq. ([1]), in which the 
left-hand side is replaced with 0, by using the imaginary- 
time method [i^. The wavefunctions are decomposed in 
a cartesian basis with a mesh-size unit A x = 0.8 fm [io^. 
The encapsulating box has to be large enough so that the 
tails of (pi are not significantly affected by the hard-box 
boundary condition. The HF calculations are converged 
for the Cf nucleus with 16 steps of A x from the center 
of the nucleus. 

The dynamical calculations for central collisions are 
performed in a half-box with = 96, Ny = 32, and 
Nz = 16 mesh-points along the x, and z-axis, respec- 
tively. The z = plane is assumed to be a plane of 
symmetry to speed up the calculations. For non-central 
collisions, Ny is doubled to allow full re-separation of the 
fragments without spurious refiections at the box bound- 
aries. The nuclei start initially along the x axis at a 
distance I^o = 51.2 fm. Their initial velocity vectors 
are determined assuming a Rutherford trajectory and 
they are given a boost by applying a translation in mo- 
mentum space [45]. Equations (j4]) are then solved it- 
eratively using a real time propagation algorithm that 
ensures energy conservation [35j, |46|] (see also Ref. [38] 
for more details). The tdhf3d code [ioj is used with 
time step of 1.5 x 10 s for a maximum simulation 
time of 6 xlO~^^s, sufficient for contact and subse- 
quent re-separation of the fragments. Figure [1] shows the 
half-box encapsulating an example of isodensity obtained 
after contact of the nuclei in a central collision. 




FIG. 1: (color online). A 12.8 fm x 25.6 fm x 76.8 fm 
half-box used for 232rp]^^250Q£ central collisions. The 
surface represents an example of isodensity, at half the 
saturation density po/2 = 0.08 fm~^, of the fragments 
moving apart after contact. 



III. MULTINUCLEON TRANSFER IN 
232r,.jj_^25o^p AT LOW-ENERGY 

Let us now study the multinucleon transfer mecha- 
nism in the 232rpj^^250Q£ reaction at energies between 
626.6 MeV (no contact) and 1205 MeV with the tdhfSd 
code. 
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A. Definition of the relative orientations 





FIG. 2: (color online). Nucleon density in the z = 
plane at various times for a central collision of a ^^^Cf 
(initially on the left) with a '^^'^Th (right) nucleus at a 
centre of mass energy Ec.m. = 1012.2 MeV. Snapshots 
are shown from t = 2Ax 10~^^ s to 3.54 x 10~^^ s in 
time steps of 3 x 10~^^ s from top to bottom. From left 
to right, the columns represent the XX, YX and YY 
relative orientations (see text). Dark blue denotes 
densities below 0.1 fm~^ and dark red marks those 
above 0.16 fm~^. 

Both nuclei exhibit a strong prolate deformation in 
their ground state and can, in principle, take all possible 
orientation in the entrance channel. Five different rel- 
ative orientations between the nuclei, labelled XX, XY, 



YX, YY, and YZ, have been selected to study their role 
on the reaction mechanism (see top of Fig. [2] and Fig.[3a|). 
We define them according to how the elongation axes are 
angled with the collision axis (i.e., the x axis). For in- 
stance, the XX (YY) orientation involves the two nuclei 
colliding on their tips (sides), as shown by the left (right) 
column of Fig.[2j In the XY and YX orientations, contact 
occurs first between the tip of one nucleus and the side 
of the other. The first letter corresponds to the heavier 
nucleus. Thus, the central column of Fig. [2] displays the 
YX orientation, with the elongation axis of ^so^f (232 r^j^^ 
perpendicular (parallel) to the collision axis. Figure [2] 
clearly shows the importance of the initial orientation 
on the reaction mechanism. For instance, in the last 
snapshot, the fragments in the XX configuration are well 
separated while a neck is still present in the other orien- 
tations, showing a shorter contact time in the XX case. 
The internal density and shape evolutions also depends 
on the orientation, going from strong fiuctuations in the 
XX orientation to a smooth evolution in the YY one. Fi- 
nally, the YX orientation produces the heaviest element 
(left fragment), corresponding to a transfermium nucleus. 



B. Multinucleon transfer in central collisions 

The process of standard quasifission is usually domi- 
nant in reactions with heavy nuclei where nucleons are 
transferred from the heavier to the lighter nucleus. As 
the dinuclear system is electrostatically unstable, it then 
separates into two fragments, with an increased mass 
symmetry. The production of transfermium nuclei in 
232rpj^^250Q£ implies that a product nucleus has to have 
more mass than either of the original two. Quasifission 
must either act in reverse, with nucleon transfer from 
the lighter to the heavier nucleus, called inverse quasi- 
fission (IQ), or overshoots so that '^^'^Th attains enough 
nucleons to end up heavier than ^^^Cf, which we define 
as swap-IQ. 

To count the number of protons, Z/., and neutrons, 
Nf. , in the fragment i (i = 1 , 2 as all the present calcu- 
lations show only two fragments in the exit channel), an 
integration of the corresponding densities is performed in 
space regions where the total density exceeds 0.001 fm~^ 
at the last iteration time. Applying this procedure at 
the initial time step excludes approximately 0.7 bound 
neutrons per fragment and no proton. The remnant neu- 
trons are found in the tails of the wave functions. Thus 
this value is added to the integration of density to obtain 
a correct estimate of Nf. . Note that the TDHF evolution 
is stopped before the fragments reach the walls of the box 
to avoid nucleon emission due to unrealistic rebounds. 

The dependence on beam energy of the proton and 
neutron numbers in the ^^^Cf-like fragment are plotted 
in Figs. [3b] and [3cl respectively, for each initial relative 
orientation. The dashed lines represent the change in nu- 
cleon value required to end either as a ^so^f q^. 232 rpj^ 
in the exit channel. Events lying between these two 
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comes into contact with the side of 250q£ absorbed 
(see the sketch of the exit channel in Figl3a|). Note that 
the rapid increase of the number of transfered nucleons 
around Ec.m. — 1200 MeV can be attributed to strong 
dynamical fluctuations of the internal density, modifying 
the breaking point of the dinuclear system (see Ref. (31j 
and section UlI Pp , rather than standard transfer where 
the flux of nucleons occurs with a smooth change of the 
shape of the fragments. 

Focusing on the low-energy range {Ec.m. < 1000 MeV) 
of the YX orientation in Fig.[3l the most massive nucleus, 
corresponding to ^^^Lr, is formed at 915.8 MeV. This 
corresponds to three neutrons heavier than the heaviest 
Lawrencium isotope found experimentally to date p^of . 
Furthermore, this corresponds only to the expected cen- 
ter of the fragment mass and charge distributions for 
this particular orientation. Taking into account par- 
ticle number fluctuations in the fragment (which are 
known to be underestimated in TDHF [47]) would lead to 
more neutron-rich nuclei in the tail of the fragment mass 
distribution. Finaly, let us note that heavy fragments 
are also produced by swap-IQ in the XY orientation at 



950 MeV. 
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C. Fast neutron evaporation 
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FIG. 3: (color online), (a) Representation of the initial 
relative orientations and sketch of the transfer from the 

tip to the side in the XY and YX configurations. 
Variation of (b) proton and (c) neutron numbers in the 
^^^Cf-like fragment after collision with a '^^'^Th nucleus 

at varying centre of mass energies for five relative 
orientations. The dashed lines at represent ^^^Cf and 

the lines at AZ = -8 and AN = -10 mark ^^^Th. 



lines correspond to standard quasifission, while events 
above the upper and below the lower Hue are associated 
with IQ and swap-IQ processes, respectively. While most 
events are located around or between these lines, the YX 
configuration leads clearly to a strong IQ for centre of 
mass energies Ec.m. > 800 MeV. Here, the tip of ^^^Th 
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FIG. 4: (color online). The total number of neutrons 
evaporated from the 250q£ 232 r^j^ fragments for 
varying centre of mass collision energies and five 
different relative orientations. 



During the collision, a significant part of the relative ki- 
netic energy of the nuclei is transformed into internal ex- 
citation of the dinuclear system and its fragments. Thus, 
the system can emit particles, in particular neutrons, be- 
fore and after its separation [48]. On one hand, this neu- 
tron emission reduces the chance to produce neutron-rich 
nuclei, but, in the other hand, it is a cooling mechanism 
which increases the survival probability of the fragments 
against secondary fission. 

In principle, nucleon emission is a one-body process 
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accounted for in TDHF. However, the finite time of the 
calculation allows only to estimate the total number of 
emitted nucleons soon after the re-separation of the frag- 
ments (typically about 10~^^ s after the neck breaks). To 
get a quantitative insight into neutron emission, the to- 
tal number of neutrons lost to the fragments is computed 
at the end of the calculation. As the TDHF evolution 
is unitary and conserves the total number of neutrons, 
Nfot. = 294, the number of evaporated neutrons is de- 
termined from the relation N^yap. = ^tot. — ^fi — Nf^. 
Figure |4] gives the evolution of N^vap. foi" various orien- 
tations and energies. A global linear increase of emitted 
neutrons is observed with energy. In contrast, the cal- 
culations shows that no proton has been lost to either 
fragment, due to the Coulomb barrier at the surface of 
the nuclei. 

The way the number of nucleons in the fragments is 
defined in section HH BI implies that the variations of the 
neutron number in the ^^^Cf-like fragment in Fig. [3c] al- 
ready takes into account this neutron emission. For in- 
stance, in the YX case at Ec.m. = 915.8 MeV, leading 
to the ^^^Lr nucleus, approximately three neutrons have 
been evaporated, i.e., ~ 1 — 2 neutron per fragment, car- 
rying away some of their excitation energy. The subse- 
quent decay occurs by neutron and gamma emission, and 
by secondary fission. The question of the remaining ex- 
citation energy of the fragments is essential to determine 
their survival probability against fission in one hand, and, 
in the other hand, to predict which isotopes are finally 
produced. For instance, for all IQ events in the YX con- 
figuration, the final total kinetic energy of the fragments 
predicted by TDHF is of the order of - 650 MeV. Then, 
to enhance the survival probability of neutron-rich heavy 
nuclei, it may be preferable to consider a lower energy 
than the optimum one deduced by Fig. [3] [49i] . 



D. Collision time and saturation in the neck 

The multinucleon transfer is expected to be affected by 
the life-time of the dinuclear system, i.e., the time dur- 
ing which the two fragments are in contact. The collision 
time is also an important input for example in the cal- 
culations of electron-positron pair production from the 
quantum electrodynamics (QED) vacuum decay (5Ql-[52j. 
As for multinucleon transfer, the collision time between 
actinides has been recently investigated in various mod- 
els [H, m, Hi-lIl [H as well as experimentally [Isj. In 
particular, it has been shown that the collision time de- 
pends on the initial orientation [31]. Indeed, as one can 
see in Fig. [2l it is much smaller for the XX orientation. 

Following Ref. [31], we define the collision (or "con- 
tact") time as the time during which the fragments are 
in contact with a neck density exceeding one tenth of 
the saturation density, i.e., Pneck > Po/10 = 0.016 fm~^. 
Figure [5] presents the evolution of this time with energy. 
The same behavior as in the 238^^238]^ ^^^^ observed 
in the present energy range (see Fig. 2b of Ref. [31]), 
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FIG. 5: (color online). Time during which the 
fragments are in contact, for five relative orientations, 
as function of the center of mass energy. 



i.e., an increase with energy up to 3 — 4 x 10~^^ s at 
Ec.m. = 1200 MeV for all orientations except the XX one 
which exhibits a plateau at ~ 2 x 10~^^ s. Comparing 
Figs. [3] and \5\ it is interesting to note that the absolute 
value of the number of transfered nucleons and the con- 
tact times have very different behaviors. This may be 
attributed to the dynamics of the dinuclear system, in 
particular in term of its complex shape evolution (see, 
e.g.. Fig. [2]). Note that the decrease of the collision time 
at higher energy observed in [31] is outside the energy 
range of the present calculations. 

As in the uranium-uranium case, the saturation of the 
collision time in the XX orientation can be interpreted as 
an effect of the overcoming of the saturation density in 
the neck, inducing a strong repulsion between the frag- 
ments. To get a deeper insight into this effect, let us 
study the maximum density for two overlapping nuclei. 
The criterion to define that the nuclei overlap is that 
the minimum density in the neck region on the collision 
axis has to be greater than 0.14 fm~^. Fig. [6] shows the 
maximum density along the collision axis during the over- 
lap. The maximum density increases with energy above 
an energy threshold which depends on the initial orien- 
tation of the nuclei. A nucleus with an orientation X, 
i.e., with its deformation axis along the collision axis, 
overlaps with its collision partner sooner, and at lower 
energy than for the other orientations. This is why the 
energy threshold Eth. above which the maximum density 
increases is lower for XX {E^^^ 870 MeV) than for YX 
and XY (^/^^ ' - 960 MeV), which, in turn, are also 
lower than for YY and YZ {E^j^^^^ - 1010 MeV). As 



a consequence, the maximum density exceeds the satu- 
ration density in the XX at lower energy (typically for 
E > 1000 MeV), than in the other orientations. It is 
interesting to note, however, that the plateau in the con- 
tact time in the XX orientation (see Fig. [5]) starts at lower 
energy (at about ~ 770 MeV) than E^^ and that other 
dynamical effects may also play a role in the saturation 
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FIG. 6: (color online). Maximum density along the 
collision axis between 250q£ ^^^Th, for varying 
centre of mass energies and five relative orientations. 
The dashed line represents the saturation density 
po = 0.16 fm~^ and points are displayed only if the two 
nuclei overlap, as decided by a minimum central axis 
density greater than 0.14 fm~^. 

We also note that the sudden increase of IQ nucleon 
transfer for orientation YX at Ec.m. - 1200 MeV in Fig. [3] 
coincides with densities surpassing saturation in Fig. [6l 
A close look at the internal density evolution shows that 
this overcoming of po generates fast dynamical fluctua- 
tions. As a consequence, the system can break at differ- 
ent positions, inducing variations of the number of nucle- 
ons in the fragments that are not totally due to standard 
multinucleon transfer through the neck. It is important 
to note that this effect occurs only in violent collisions, 
and that the excited fragments should have a very small 
chance to survive against subsequent fission. 



E. Role of impact parameter 

We finally investigate how multinucleon transfer 
evolves with impact parameter for the YX orientation at 
Ec.m. = 915.8 MeV, where the heaviest nucleus (^^^Lr) 
is formed. As described in section ITiBi these non-central 
collisions are performed in a twice bigger box to avoid any 
spurious effect of the box before the full re-separation of 
the fragments. Like in the central collision case, the nu- 
clei start initially on the x axis. However, their initial 
Rutherford trajectory is determined for a finite impact 
parameter b. 

Figs. [Til and [7b] display the post-collision number of 
protons and neutrons in the heavier fragment, respec- 
tively. The global effect of increasing the impact param- 
eter is to reduce the number of nucleons transferred via 
IQ. This can be interpreted in terms of a reduction of the 
contact time because of the centrifugal potential. Note 
that the transfermium production is predicted to be dom- 



inant for this configuration up to 6 ^ 3 fm, corresponding 
to an angular momentum of ^ llSh. 

It is interesting to note that IQ disappears for impact 
parameters above b 4 fm. In particular, the heavy 
fragments loses about one proton, but no neutron, for 
4 < b < 8 fm. This can be understood in terms of 
charge equilibration as ^so^f slightly more proton-rich 
than 232 rpj^ ^ transfer of one proton is enough to 
equilibrate this asymmetry. In addition, with the initial 
condition for a YX orientation, a non-zero impact pa- 
rameter shifts the system at contact, going away from 
the configuration where the tip of one nucleus collides 
with the side of the other, which we identified as the 
most favorable in terms of heavy elements production in 
section IIIIBl Finally, at 6 > 8 fm, the overlap is not 
sufficient to allow any transfer of nucleons. 
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FIG. 7: (color online). Number of (a) protons Z and 
(b) neutrons N in the heavier fragment as function of 
the impact parameter at a center of mass energy 
Ec.m. = 915.8 MeV in the YX relative orientation. 
^^^Cf is represented by the dashed lines at Z = 98 and 
N = 152. 
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IV. CONCLUSIONS 

The time-dependent Hartree-Fock theory has been 
used to study the reaction mechanisms in the 
232rp]^_^250(-;£ reaction. The role of the deformation and 
relative orientation has been investigated and number of 
transfered nucleons, collision time, density in the overlap 
region and fast neutron emission have been analyzed. 

A new process of inverse quasifission has been iden- 
tified when the tip of the lighter nucleus collides with 
the side of the heavier one. In this case, nucleons are 
transfered to the heavier nucleus and new neutron-rich 
transfermium nuclei can be produced. With the present 
reaction, ^^^Lr, which has three more neutrons than the 
most neutron-rich observed lawrencium isotope, could be 
produced in this process. In addition, fluctuations in the 



fragment neutron distribution should produce even more 
neutron-rich nuclei. 
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